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In this paper we study the properties of two-qubit gates. We review the most common parameter- 
izations for the local equivalence classes of two-qubit gates and the connections between them. We 
then introduce a new discrete local invariant, namely the number of local degrees of freedom that 
a gate can bind. The value of this invariant is calculated analytically for all the local equivalence 
classes of two-qubit gates. We find that almost all two-qubit gates can bind the full six local degrees 
of freedom and are in this sense more effective than the controlled-NOT gate which only can bind 
four local degrees of freedom. 

PACS numbers: 03.67.Lx 

Keywords: quantum computation, local invariants 



INTRODUCTION 



Quantum computation is a novel information processing method in which classical information is encoded into a 
quantum- mechanical system pj, called the quantum register. In most quantum computers the quantum register is a 
■ collection of two-level systems, termed qubits. The computation is performed by the unitary temporal evolution of 
' the register, followed by a measurement. In order to execute a quantum algorithm, one has to be able to generate the 
. required unitary propagators that are usually referred to as quantum gates. 

o 



It has been shown that almost any fixed two-qubit gate together with arbitrary single-qubit gates is universal 0,01, i-e., 
any n-qubit gate may be constructed using only a finite number of these gates. Conventionally, the elementary gate 
library is chosen to consist of the single-qubit rotations R x ,R y ,R z and the controlled-NOT gate (CNOT). However, 
in many realizations, the CNOT is not the natural choice for the entangling two-qubit gate. Recently, an optimal 
construction of an arbitrary two-qubit gate using three CNOTs and 15 single-qubit rotations has been introduced 0|. 
In addition, constructions for the double-CNOT (DCNOT) 0, the controlled-unitary gates |g] and the so-called súper 
q-i controlled gates Q have been published. A construction using only two applications of the B gate has been introduced 
in Ref. 0, and in Ref. it is shown that no other construction using only two applications of a fixed two-qubit 
. £h ! gate exists. Extensions to the n-qubit case are mainly uninvestigated. However, several CNOT-based constructions 
with 0(4") asymptotic behaviour exist, the best of which have CNOT counts of twice the highest known lower 

bound Ü|. 

_ _ i 

In many of the proposed realizations for quantum computers the individual qubits are fully controllable, whereas the 
interqubit interactions are often fixed. In addition, single-qubit operations tend to be considerably faster to implement 
than multiqubit operations. This is why it often makes sense to study the local equivalence classes of multiqubit gates 
instead of the gates themselves. Two gates are considered equivalent if they can be converted to each other using 
only local operations, i.e., tensor products of single-qubit gates. The equivalence classes are characterized by local 
invariants, which are quantities that are not affected by local operations. 

In this paper we briefly review the currently used parameterizations for the local equivalence classes of two-qubit 
gates and point out their equivalence. We then introduce a new discrete local invariant which describes the number of 
local degrees of freedom a gate can bind. Finally, we calculate the value of this invariant for all the local equivalence 
classes of two-qubit gates. 
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II. LOCAL EQUIVALENCE CLASSES OF TWO-QUBIT GATES 



An n-qubit quantum gate k is said to be local iff it consists solely of single-qubit rotations: k £ SU(2)® n —: L„. Two 
n-qubit gates U\,Ui € SU(2 n ) are said to be locally equivalent iff U2 = kiUik 2 , where ki, k 2 € L n . This constitutes 
an equivalence relation, which we denote by Ui ~ C/2. 

Using the theory of Lie groups it can be shown flUÍÏ^ that any two-qubit gate U € SU (4) can be decomposed using 
the Cartan decomposition as 



U = kiAk 2 = fei exp ( -(ci<t x 



(1) 



where o", denote the Pauli matrices, ki, fe € L2 and ci, C2, C3 € R. The matrix A is a member of the Cartan subgroup 
of the decomposition and carries all the nonlocal properties of the gate U. Hence the local equivalence classes of two- 
qubit gates can be parameterized by the three scalars [ci,c 2 ,C3], known as canonical parameters. This is a minimal 
set of parameters since the group 5*1/(4) is 15-dimensional and the local rotations eliminate 2 x dim(S'C/(2)® 2 ) = 12 
degrees of freedom thereof. The canonical parameterization is visualized in Fig.Q] The tetrahedron OA1A2A3 in the 
figure is called a Weyl chamber. It is defined by the inequalities ir > c\ > c 2 > C3 > 0, 7r — c\ > C2. The Weyl chamber 
contains all the local equivalence classes of two-qubit gates exactly once, excepting the fact that the triangles LAxA 2 
and LOA2 are equivalent. 




A,= [n/2, n/2, n/2] 



[n/2, n/2, 0] 



Figure 1: Weyl chamber. Points O and Ai correspond to the identity gate /, A3 to the SWAP gate, L to the controlled-NOT 
gate (CNOT) and A 2 to the double controlled-NOT gate (DCNOT) [ÏJ. 



The matrix 




(2) 



is the transformation from the Standard basis of states {|00) , 1 01) , 1 10) , |11)} into the Bell basis, also known as the 
màgic basis We use the lower index B to denote the change of basis: Ub '■= QUQ^. The màgic basis has 

the special property that local gates expressed in it are orthogonal. In other words, conjugation by Q is a group 
isomorphism between SU{2) ® SU{2) and SO{4). Furthermore, it renders our chosen Cartan subgroup (generated by 
o~ x <8> &x, cf v ® <j y and a z ® a z ) diagonal. These two properties enable us to calculate the canonical parameters of any 
given SU(A) gate U — k\Ak2- The parameters are obtained from the spectrum of the matrix M(U) := U^Ub which 
is given by 



A (M(U)) = | 



— J p i{ci+c 2 -c 3 ) i(ci-c 2 +c 3 ) i(-ci+c 2 +c 3 ) -i(c 1 +c 2 +c 3 ) 



(3) 
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Table I: Vàlues of the canonical and Makhlin invariants for some common gates. SPE denotes a special perfect entangler [3.lïfil|. 



Ref. presents an algorithm for extracting the canomcal parameters Cj from this spectrum in a convenient way 
although it uses a slightly different notation. The equivalence of the methods becomes apparent using the equality 
Q T Q = —a y ® (j y , since 

A (M(U)) = X ((QUQÏfQUQÏ) = X ((<r y ® a y )^U T (a y ® a y )U) = X (U(a v ® cr y )U T (a y <g> cr,,)) = A(C/f>). (4) 

Ref. Ü[ presents another system of invariants, namely the characteristic polynomials x[72(£0L where 72(1/) = U(a y ® 
a y)^ T ( a v ® fj/)- They are completely equivalent to the canonical parameters since the characteristic polynomial 
x[72(C/)] carries exactly the same information as À (M(Z7)) = A(72(C/)). 

Another useful parameterization for the two-qubit local equivalence classes is provided by the Makhlin invariants G\ 
and G2 For a gate E7 S Í7(4), they are defined as 

_Tr 2 M([/) _ Tr 2 Af(C/)-TrM([/) 2 

Gl " 16detC7 ' G2 ~ 4deTÏ7 ' (5) 

The Makhlin invariants are by far the easiest ones to calculate. They, too, provide the same information as the 
previous invariants since A (M(U)) is fully determined by them. G\ may be complex but G2 is always a real number, 
which leads to three real-valued invariants. If U is represented as in Eq. {JJ, the Makhlin invariants reduce to ^2] 

gi := Re G\ = cos 2 c\ cos 2 C2 cos 2 C3 — sin 2 c\ sin 2 C2 sin 2 C3, 

,92 := Im G\ = — sin 2c\ sin 2c2 sin 2c3, (6) 

Q O O O O 

53 := G2 = 4 cos ci cos C2 cos C3 — 4 sin ci sin C2 sin C3 — cos 2ci cos 2c2 cos 2C3. 

Example vàlues of the invariants of different gates are given in Table Q| The set of all the two-qubit gate equivalence 
classes in the Makhlin parameter space is presented in Fig. |2j The surface is given by the equations 

gi — cos 2 s cos 4 1 — sin 2 s sin 4 1 

g 2 = isin(2s)sin 2 (2í) (7) 
53 = 4.gi - cos(2s) cos 2 (2í), 

where s € [0, ir], t S [0, n/2]. The surface and the inside of the object correspond to the surface and the inside of the 
Weyl chamber, respectively. 
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Figure 2: Weyl chamber in the coordinates of the Makhlin invariants. 



where j runs over the 3n local generators of SU(2 n ), to denote a fc-parameter family of n-qubit local gates. It is 
defined by the function a : M fe — > R 3n . The generators Xj are normalized such that they are orthonormal with respect 
to the inner product (X, Y) := Tr (X^Y). 

A gate U 6 SU(2 n ) is said to leak k local degrees of freedom iff there exist nondegenerate functions a and b such 
that 



UL k n ( a ,6)=L k n (h 1 0)U V6e 



(9) 



A gate binds the local degrees of freedom that it does not leak. We define a function rj : SU(2 n ) — > N to indicate the 
number of local degrees of freedom that an n-qubit gate U binds. We always have maxíy < 3n, i.e., at most three 
degrees of freedom for each qubit. 

Assume now that the functions a and b satisfy Eq. ^ for the gate U. For a gate V = k\Uki, where fei, ki € L n , we 
obtain 



V 



klL k n (a,e)k 2 ] = k 1 Uk 2 klL k n (aL,e)k 2 = fciL*(b,^ Uk 2 = \kxL k n (h,e)k\ 



V. 



We also have 



3rt 



3n \ 

hL k n (b,6)kt = exp \Y,bj(9) Ad^Xj ] = exp i^W)^ = L k (h,9), 



(10) 



(11) 



since Ad(g) is a linear bijection and ki is a local gate. If b is nondegenerate then so is b. A similar argument naturally 
holds for k 2 L k (a.,9)k 2 , which yields VL k (à,6)L k (h,9)V and proves that 77 is indeed a local invariant. 

Equation Q is equivalent to 



3n 



C/exp ( J2 a i(9)Xj C/ f = exp l^ajfyMfflXj = exp ^ò fe (0)A fc 
u'=i / \i=i / Vfe=i 



(12) 
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This is fulfilled if 

3n 3n 

J2 *M MU)x 3 = J2 h(0)x k . (13) 

j-\ k=l 

Now, as we take a sidewise inner product (•, Xi) with each of the 4™ — 1 generators of SU (2™), we obtain equivalently 

3n 3n 

Y,W ijaj (&)=Y,b k (6)S ki , i = l,2,...,4"-l, (14) 
i=i fc=i 

where = Tr (t/XÍí/tx^ . The generators Xj are antihermitian and U is unitary. This implies that the elements 
Wij are real. Written in matrix form this is 

^=Wa(í)=^j VÍeR* (15) 

where S ]R 3nx3n i Wn € J^ 4 "- 1 - 3 ™)* 3 ™ an( j ^ e mc [i ces £ and iV stand for local and nonlocal, respectively. Hence, 
we must have a(9) S ker Wn for all vàlues of 9. Moreover, since b must have the same dimensionality as a, the 
component of ker Wn parallel to ker Wl must be disregarded. Using the rank-nullity theorem we finally obtain 

V (U) = 3n - dim(ker W N ) + dim(ker W L n ker W N ). (16) 



IV. r/ FOR TWO-QUBIT GATES 



For the set of two-qubit gates U € SU (4), maxtj < 6. It is obvious that 77(1) = and 77(SWAP) = since all local gates 
and hence all local degrees of freedom may be commuted through these gates. It is also known that ?y(CNOT) = 4 
and ry(DCNOT) = 4. The result for CNOT is obtained by combining the commutation properties of CNOT with 
the Euler rotations R z and R x and the fact that an arbitrary two-qubit gate may be implemented using at most 
three CNOTs i E El El- 

Similar arguments for the DCNOT are presented in Ref. |5j], includingthe explicit 
implementation of an arbitrary two-qubit gate using three DCNOTs. Also, from the construction of Ref. pj, it is clear 
that rj(B) > 5. Apart from such observations, no explicit calculations for 77 have been presented in the literature so 
far. 



We will now proceed to derive an analytical expression for 77 for an arbitrary two-qubit gate. Because r\ is a local 
invariant, it is enough to consider gates of the type 

/ C 3 





A = exp ( -{ciacc 
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(17) 



which represent all the nonlocal equivalence classes. The calculation of the elements of Wl and Wn is straightforward. 
Calculating the matrix exponential and simplifying the expression using elementary trigonometric identities results 
in 
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[ci,C2,C3] Set in the Weyl chamber r\ 

[0, 0, 0] = [tt, 0, 0] O, Ax 

[tt/2, tt/2, tt/2] ^ 3 

[x,a;,x],i / 0,x £ tt/2 OA 3 \ {O, A 3 } 3 

[tt - x,s,a;],s / 0,a; / vr/2 A1A3 \ {Ai, A 3 } 3 

[k,0,0] = [tt-s,0,0],x/0 OA 1 \{0,A 1 } 4 

[n/2,ir/2,x],x^ir/2 A 2 A 3 \ {A 3 } 4 

[x,x,y],x=£y,x=£ix/2 OAiA 3 \ {OA 3 , AíAs} 5 

[n-x,x,y],x=£y,x£ir/2 AiA 2 A 3 \ {AiA 3 , A 2 A 3 ] 5 

{All other points} {All other points} 6 



Table II: 17, or the number of local degrees of freedom bound, for the local equivalence classes of two-qubit gates. 



where the non-zero elements are 

1 M11 ^1,2) _ [COSC2COSC3 sinc2sinc3\ .-1 _ | n !,i n i.2 \ _ ( sin ci cos C3 — cos C2 sin C3 



J^l "1,1 "1,2 ™" u**»^»™*^ 

'21 '22/ \sinc2sinc3 COSC2COSC3/ ' \ n 2 1 n 2 2 J \ — cosc 2 sinc3 sin C2 COS C3 



^2 Mi,i '1,2 ] _ / cos Ci cos C3 sincisinc 3 \ ^ 2 ( n ï,i n i.2j _ / — sinciCosc3 COSC1SÍ11C3 

V2 1 '22/ \sinc1sinc3 cos ci cos C3 I \ n 2 1 "22/ \ cosciSÍnc3 — sin c\ cos C3 

^ 3 / i if i2 ] _ /cos ci cos C2 sincisinc 2 \ ^3 [ n i,i n i.2j _ ( sinciCOSc 2 — coscismc 2 

V2 1 ^2.2/ \SÍnClSÍnC2 COS Cl COS C2 I \ n 21 n 2.2l \ — COS Cl SÍn C2 SÍll Cl COS C2 



(19) 



From Eqs. Ill8fl - I|19|l it is seen that Eq. lllót decomposes into six separate equations: 

L*( ai ) = ( h h ), N i ( at ) = (°X ï = 1 < 2 ' 3 - ( 2 °) 
\a»+3/ \bi+3j \ai+3j \0J 

Each block L' 1 produces a two-dimensional null space iff all the elements of U equal zero. A one-dimensional null 
space is formed iff detL 1 = cos(cj + Cfc) cos(cj — c k ) — 0, where e^fe = 1, but N 1 ^ 0. Similarly, each block 
produces a two-dimensional null space iff all the elements of N l equal zero, and a one-dimensional null space iff 
det iV' = (-l)' i+1 sm(c 3 + c k ) sm{c 3 - c k ) = 0, where e ljk = 1, but N l ^ 0. 

Taking into account the correlations among the elements of the matrices U and N l we find that in the two-qubit 
case ker Wl H ker Wn = {0} always. Thus we have rj = 6 — dim(ker Wn) and the number of local degrees of freedom 
leaked is given by the nullity of Wn- The results for all the possible vàlues of [ci, C2, C3] are collected in Table lïïl One 
notices that everywhere inside the Weyl chamber 77 reaches its maximum value of 6. At the vertices O = A\ and A3 
i] = 0, on the edges between them 77 = 3, on the edges OAi, A2A3 rj — 4 and on the faces OA\ A3, OA2A3, A1A2A3 
t] = 5. 

The number of local degrees of freedom that the gate U leaks is obtained as the number of pairs of equal eigenvalues 
Ai in the spectrum of the matrix M(U), presented in Eq. (EJ. In other words, any n-fold eigenvalue of M(U) indicates 
n(n — l)/2 local degrees of freedom that pass through the gate U. Translated to the language of the Weyl chamber, 
each Weyl symmetry plane the point [c\, C2, C3] touches causes the gate to leak one local degree of freedom. 
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V. CONCLUSION 



In this paper we have introduced a new local invariant rj for quantum gates, indicating the number of local degrees 
of freedom a gate can bind. Furthermore, we have analytically calculated the value of this invariant for all two-qubit 
gates. We have found that almost all two-qubit gat es can b ind the full six local degrees of freedom. However, most of 
the commonly occurring gates such as CNOT or VSWAP are exceptions to the rule, performing much worse in this 
sense. 



The meaning of r\ is illustrated by considering the lower bounds on gate counts for a genèric n-qubit circuit. Let 
the gate library consist of all one-qubit gates and a fixed two-qubit gate U. Then almost all n-qubit gates cannot be 
simulated with a circuit consisting of fewer than 



N„ = 



4" - 3n - 1 

— W) 



(21) 



applications of the two-qubit gate. This result is a straightforward generalization of Proposition III. 1 in Ref. The 
gates binding the full six degrees of freedom are thus expected to be the most efficient building blocks for multiqubit 
gates. 
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Appendix A: MATHEMATICAL PREREQUISITES 



The Lie àlgebra g of a linear Lie group G < GL(n, K) is the set 

g:={X e K nxn \ exp(íX) e G Vt G K}. (Al) 

In can be shown that g is a real vector space spanned by the generators of G. For example, the Lie àlgebra su(n) of 
the group SU(n) consists of all the n x n complex antihermitian traceless matrices. 

The adjoint representation of a Lie group G, Ad : G — > Aut(g), is a group homomorphism defined by 

Ad(g)X := gXg- 1 ( g eG,Xeg). (A2) 
It behaves in a rather simple way in exponentiation: 

exp {Ad(g)X) = g exp (X) g~ l for all g £ G, X e g. (A3) 
Also, if we define an inner product (X, Y) := Tr (AT^F) for g, we find that it is preserved by the adjoint representation: 



(Ad(g)X, Ad(g)Y) = (X, Y) for all g G G, X, Y e g. 



(A4) 
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As a concrete example, the adjoint representation keeps orthonormal bases of g orthonormal. 
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